Noise correlation spectroscopy of the broken order of a Mott insulating phase 
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We use a two-color lattice to break the homogeneous site occupation of an atomic Mott Insulator of 
bosonic 87 Rb. We detect the disruption of the ordered Mott domains via noise correlation analysis 
of the atomic density distribution after time-of- flight. The appearance of additional correlation 
peaks evidences the redistribution of the atoms into a strongly inhomogeneous insulating state, in 
quantitative agreement with the predictions. 
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Ultracold atoms in optical lattices are optimum candi- 
dates for quantum computation [l[ and quantum simula- 
tion [2], |3| thanks to the possibility to control and access 
the atomic state. They have demonstrated great capabil- 
ities in reproducing ideal solid state systems, a remark- 
able example being the realization of the Mott Insulator 
(MI) phase 0,0]. Furthermore, experiments are now be- 
ginning to access a rather unexplored class of quantum 
systems involving correlations and disorder [|| . The need 
for effective diagnostic techniques being able to observe 
these quantum phases and compellingly establish their 
nature has become notable. Few years ago Altman et al. 
suggested in ref. 7] the use of Hanbury Brown and Twiss 
(HBT) spatial interferometry to probe hidden order in 
strongly correlated phases of ultracold atoms. This pro- 
posal has been recently employed to experimentally in- 
vestigate quantum statistically correlated lattice systems 
[1, B fioj and momentum pair-correlated fermions after 
molecule dissociation [ll[. Meanwhile, theoretical works 
have extended the use of noise interferometry to study 
disordered systems [HI, [H| and to observe a predicted 
supersolid phase [Til ]. 

In this work we break the ordered domains of a one- 
dimensional MI in a controlled way, by means of a sec- 
ond periodic potential superimposed on the main lattice. 
The MI is a state of matter in which strong interactions 
force bosonic particles to localize at the lattice sites, in 
the homogenous case with exactly the same number of 
particles per site. This energetic stiffness of the phase 
remains even if, as it happens in usual experiments with 
cold atoms in optical lattices, a "mildly" varying trapping 
potential is superimposed. This variation of site en ergy 
produces the well known "wedding cake" structure [15| . 
in which the system rearranges regularly alternating MI 
domains with different occupation number and superfluid 
islands. Locally the main properties of the homogeneous 
system are still preserved, and, as confirmed by exper- 
iments [B|, [r| the global system has mostly an ordered 
insulating character. A more intimate alteration of the 
MI structure can be instead obtained by a non-monotone 
scrambling of the energy offsets of the sites, employing a 
two-color lattice. 

As earlier shown in ref. @, [13], the superposition of 
an additional weak lattice to a MI state is enough to 
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FIG. 1: In a MI domain the atoms are homogeneously dis- 
tributed pinned by strong interactions at the lattice wells. 
The superposition of a second weaker lattice with different 
color (shown in the lower row) breaks the translational in- 
variance leading to a redistribution of particles among lattice 
sites resulting in a final non-homogeneous density configura- 
tion. The change in the site occupation is the main effect 
of the additional lattice, as the position of the wells has not 
varied significantly. 



profoundly affect the excitation spectrum of the system. 
When this second lattice is strong enough, we expect the 
atoms to redistribute over the lattice breaking the or- 
dered structure of the MI. In this Letter we describe an 
experiment where we have detected this redistribution 
using the noise interferometry technique in a yet unex- 
ploited way. In the experiments performed until now, 
the noise correlation has only given information on the 
strongly correlated character of the MI and on the spa- 
tial ordering of the in-trap density distribution, namely 
the fact that atoms are located at equidistant positions 
in space. In this work we further demonstrate the pow- 
crfulness of this diagnostic technique by showing that it 
can be used to detect modifications in the ordered site 
filling typical of a MI state, whereas the density distribu- 
tion after time-of-flight (TOF) does not show qualitative 
differences. 

In the experiment we load a 87 Rb Bose-Einstein con- 
densate (BEC) containing typically 3 x 10 5 atoms in a 
2D optical lattice (confining lattice) with height s t = 40, 
in units of the recoil energy Em — h 2 k\/2m ~ h x 3.33 
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kHz, where 2ir/ki = 830.3(1) nm is the wavelength ol 
the lattice light we derive from a Titanium-Sapphire laser 
and m is the atomic mass. In this way we create a bidi- 
mensional array of ID atomic tubes with maximum 200 
atoms, among which the tunneling is highly suppressed 
on the timescale of the experiment. Simultaneously, we 
apply along the axis of each atomic tube a third optical 
lattice (main lattice) with the same spacing and height 
si = 16. 

The system can be described by the Bose-Hubbard 
Hamiltonian Hi 
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where &j and are respectively the annihilation and 
creation operators of one boson in the j-th site of the 
lattice and fij is the number operator. The relevant en- 
ergy scales in the process are the tunneling energy J, the 
on-site repulsive atomic interaction U and the site ener- 
gies Cj determined by the external trapping potential. In 
our realization, the system is deeply in the MI phase: the 
interaction energy is U — h x 2.7 kHz and the tunneling 
along the tubes direction (that we will call x direction 
in the following) is J = h x 20 Hz. The lattice beams 
have a gaussian shape with a typical waist of ~ 200 /im 
and cause an additional external harmonic confinement. 
The resulting trapping frequencies are u x — 2tt x 77 Hz 
along the direction of the tubes and u± = 2ir x 110 Hz 
orthogonally. 

By adding a weaker second lattice with incommensu- 
rate wavelength 2ir/k2 = 1076.9(1) nm with respect to 
the main one, along the same direction, we manipulate 
the third term in the Hamiltonian |T]) according to: 
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with s 2 the height of the second lattice expressed in units 
of the recoil energy Er 2 = (tt 2 ) 2 /(2m) ~/ix 1.98 kHz. 
This additional lattice has the main effect of producing 
a non-periodic modulation of the energy minima over a 
length scale of 1.8 /im, thus breaking the translational 
invariance of the main lattice. As far as s 2 <C Si, the pa- 
rameters J and U and the spatial position of the energy 
minima Xj arc only slightly affected by the additional lat- 
tice [l7|- Increasing sg, the ordered atomic distribution 
in the Mott domains is gradually affected and atoms are 
expected to rearrange in a non uniform pattern as shown 
in Fig. [T] 

We probe the new density distribution using the noise 
correlation analysis at fixed height of the main lattice 
si = 16 and various second lattice intensities in the range 
< s 2 < 9. We load the BEC into the final optical 
lattice configuration, switching on all the lattices at the 
same time by means of an exponential ramp with time 
constant 30 ms and duration 140 ms. After holding the 
atoms for 30 ms in the lattices, we let the atomic system 
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FIG. 2: Left column: averaged images of the atomic den- 
sity distribution after 20 ms time-of-flight with si = 16 and 
varying S2. Right column: correlation functions. For S2 < 3 
peaks at a distance d = 2hk\t/m appear reflecting the main 
lattice periodicity. When S2 > 5 we observe the growth of 
additional peaks corresponding to the periodicity of the sec- 
ond lattice (fe) and to the beating between the two standing 
waves (k± — k2)- 



expand by suddenly switching off all the confining po- 
tentials. After 20 ms of ballistic expansion, we image the 
cloud by standard destructive absorption technique. The 
probe beam has an intensity / = 0.5 mW/cm 2 ~ 0.3/ sa t 
where I sa t is the saturation intensity of the Z?2 transi- 
tion. We illuminate the atomic cloud for 50 /is and then 
we image the transmitted light on an interline CCD cam- 
era. 

From the obtained images we calculate the normalized 
density-density correlation function analogously to what 
is done in ref. [H, by: 



J dr(n(r)h(r + d)) t 
J dr(n(r)) t {n{r + d))t 



(3) 



where n(r) is the column density detected at the posi- 
tion r reached by the particles in a time t after releasing 
from the trap, d is the separation between two detected 
positions in the expanded atomic cloud and () refers to 
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a statistical averaging over a set of several images taken 
in the same conditions [l8[ . A Butterworth high-pass fil- 
ter is applied on the correlation signal to eliminate the 
gaussian background appearing as a consequence of shot- 
to-shot ±25% fluctuations in total atom number. A sig- 
nificative increase in the visibility of the correlation peaks 
has been obtained by averaging the final image with itself 
after flipping around the central vertical axis. This pro- 
cedure is justified by the symmetry of the lattice system 
under inversion x — > —x. 

Non-trivial correlations have shown to be present in 
the second-order correlation function £7(d) of the density 
distribution of a bosonic atomic cloud freely expanding 
from a MI state, with peaks emerging spaced by the dis- 
tance 2hkt/m being k the wavevector of the lattice 
laser light. Figure [5] shows images of the density distribu- 
tion after TOF (left column) and the correlation function 
(right column) both averaged on 5 -j- 6 sets of 40 -j- 80 im- 
ages for different values of s 2 . The TOF images exhibit 
an unstructured gaussian profile and no significative dif- 
ferences increasing the height of the second lattice (apart 
from an increase in size along the direction of the two- 
color lattice). When s 2 < 3 the measured correlation 
signal presents peaks clearly visible in correspondence of 
the periodicity of the main lattice (k\ peaks) and no sig- 
nificative alteration with respect to the pure MI case. At 
s 2 — 5 additional peaks are detectable at the periodicity 
of the second lattice (k 2 peaks), in the central row of the 
noise correlation signal. Similar correlation peaks have 
been predicted to appear first for hard-core bosons in a 
quasiperiodic potential and then also for soft-core 
bosons [l3| . The additional peaks strengthen, becoming 
clearly visible at s 2 = 9, where also the beating between 
the two standing waves {k\~k 2 peaks) is visible in the up- 
per and lower rows. Note, that the additional peaks only 
appear along the horizontal direction where the second 
lattice modifies the site occupation. The ki — k 2 peaks 
should be visible also in the central row, but they are 
obscured by a central black ring artifact deriving from 
the application of the high-pass filter. Because the sec- 
ond lattice does not appreciably shift the position of the 
potential wells, in the range of parameters we used, the 
appearance of the additional peaks arises from a differ- 
ent effect, namely the change in the local atom site oc- 
cupation induced by the second lattice. Atoms rearrange 
in the main lattice wells searching for the lowest energy 
configuration driven by the interplay between the second- 
color modified site energies and the interactions. 

In order to obtain a deeper insight, we have devel- 
oped a theoretical model working in the zero-tunneling 
limit which confirms this interpretation. In this limit 
the Hamiltonian in equation (TjQ) is diagonal on the ten- 
sor product of number Fock states in each lattice site, 
describing a MI state. The ground state is given by the 
distribution of atoms minimizing the total energy of the 
system. The lowest energy configuration is computed, 
fixing the number of atoms, and subsequently noise cor- 
relations are calculated according to ref. Q. In Fig. [3] 
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FIG. 3: The ground state of the system is calculated by diago- 
nalizing Eq. [T]m the limit J — > 0. (a) Calculated atom den- 
sity distribution at si = 16 and S2 = 0, 2, 5, 9. (b) Noise cor- 
relation signal computed for each of the ground states given 
in (a), (c) Horizontal cross section of the peaks in the noise 
correlation signal. 



we report the result of this calculation for different val- 
ues of the second lattice intensity s 2 = 0,2,5,9. The 
picture shows in column: (a) the density distribution for 
the central atomic tube, (b) the 2D correlation signal and 
(c) its cross section along the horizontal line. At s 2 = 
the system shows the typical wedding cake structure, due 
to the presence of the harmonic trapping potential, with 
MI domains of N — 1,2,3,4 atoms per site. The cor- 
responding noise signal exhibits the peaks at periodicity 
k\ as expected. When s 2 = 2 the MI order begins to be 
destroyed starting from the edge of each domain, where 
the site energies scrambling is more efficient due to the 
different atom occupation of next-neighboring sites with 
almost the same energy. The noise signal shows weak 
peaks emerging at the periodicity k 2 and k\ — k 2 , as can 
be seen in the cross section picture. Increasing the height 
of the second lattice to s 2 — 5, the density profile gets 
more and more inhomogeneous, affecting also the bulk 
of the domains. The extra peaks are now clearly visi- 
ble in the noise signal. At the maximum value s 2 = 9, 



4 



0.8 



■B 0.6 

03 



jB 0.4 

CO 
CD 
Cl 



O 



r 0.2 



0.0 



* 



2 4 6 8 

second-color lattice height s 



10 



FIG. 4: The ratio between the height of the hi peak and the 
height of the ki peak is reported as a function of sa, while si 
is fixed at the value of 16. Each point is an average of 5 -r- 9 
values from different sets of images. The grey curve derives 
from calculation by a zero-tunneling model for a total atom 
number ranging from 2 x 10 5 to 4 x 10 5 . 



sites with N — 5 occupation appear together with empty 
sites even at the trap center and the wedding cake struc- 
ture is completely lost being washed away in favor of the 
second-color driven density distribution. The cross sec- 
tion graph shows the appearance of higher-order beat- 
ing peaks. Note that the peaks at &2 are slightly higher 
than the ones at k± — k%. We have verified in our nu- 
merical simulations that this is the only effect due to the 
small modification of the minima position induced by the 
strong second lattice (the relative deviation with respect 
to the main lattice spacing is ~ 4% at s 2 =9). 

On a qualitative level, the theoretically derived signals 
show the same behavior as the experimental observations, 
see Figj2] For a more quantitative analysis we compare 
in FigfJ]the theoretical and the experimental results for 



the ratio between the height of the hi peak and the main 
ki peak for different values of s 2 . Each point in the graph 
results from the weighted average of 5 — 9 experimental 
values. The dotted curve shows the theoretical predic- 
tion for N = 3 x 10 5 atoms. The grey zone takes into 
account possible systematic errors in the measurement of 
the absolute number of atoms, showing the results of the 
calculations for a number of atoms ranging from 2 x 10 to 
4 x 10 5 . We find good agreement between the experimen- 
tal data and the results of this zero-tunneling model with 
no free parameters. Results compatible with our exper- 
imental findings have also been obtained by T. Roscilde 
[l3T ] with a Quantum Monte Carlo calculation performed 
for a one dimensional system with N = 100 atoms, cor- 
responding to one atomic tube of our realization. 

In conclusion, we have destroyed in a controlled way 
the Mott Insulator ordered structure by using an addi- 
tional lattice at a different incommensurate wavelength 
and we have monitored its progressive degradation via 
quantum noise interferometry. We have demonstrated 
that the noise correlation technique can be employed 
to detect modifications in the occupation of the lattice 
sites. Furthermore, a zero-tunneling model shows good 
agreement with the experimental data. Future perspec- 
tives of extending our theoretical model to keep into ac- 
count finite tunneling appear interesting, especially for 
the search of possible signals that could give information 
on characterizing complex quantum phases like the Bose 
Glass d, [H, HI 
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